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Abstract 

We analyze the dynamics of N interacting spins (quantum register) collectively 
coupled to a thermal environment. Each spin experiences the same environment 
interaction, consisting of an energy conserving and an energy exchange part. 

We find the decay rates of the reduced density matrix elements in the energy 
basis. We show that if the spins do not interact among each other, then the fastest 
decay rates of off-diagonal matrix elements induced by the energy conserving inter- 
action is of order N'^ , while that one induced by the energy exchange interaction 
is of the order N only. Moreover, the diagonal matrix elements approach their 
limiting values at a rate independent of TV. For a general spin system the decay 
rates depend in a rather complicated (but explicit) way on the size N and the 
interaction between the spins. 

Our method is based on a dynamical quantum resonance theory valid for small, 
fixed values of the couplings. We do not make Markov-, Born- or weak coupling 
(van Hove) approximations. 



1 Introduction 

Description of the problem. We consider a qubit register of size N whose Hamil- 
tonian is of the form 

N N 

H^=Y,-h,StS^^+Y,B,Sl (1.1) 

i,j=l j=l 

where the Jij are pair interaction constants that can take positive or negative values, 
and Bj > is an effective magnetic field at the location of spin j {Bj = ^jB], where 
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h is the Planck constant, 7 is the value of the electron gyromagnetic ratio and Bj is 
an inhomogeneous magnetic field, oriented in the positive z direction). Also, 
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(1.2) 



is the Pauli spin 1/2 operator; Sj is the matrix acting nontrivially only on the j-th 
spin. The environment R is modelled by a bosonic thermal reservoir whose Hamiltonian 
is 

Hr= f a*{k)\k\a{k)d^k, (1.3) 

where a*{k) and a{k) are the usual bosonic creation and annihilation operators sat- 
isfying the canonical commutation relations [a{k),a*{l)] = 5{k — I). It is understood 
that we consider R in the thermodynamic limit of infinite volume, fixed temperature 
T = 1//? > 0, in a phase without Bosc-Einstein condensate. 

We consider a collective coupling: the distance between the N qubits is smaller 
than the correlation length of the reservoir and consequently each qubit feels the same 
interaction with the latter. The collective interaction between S and R is given by the 
operator 

N N 

V = Xivi + X2V2 = Ai ^ 5| ® </)(5i) + ® (t>{g2). (1.4) 

j=i j=i 
Here, <p{g) is the field operator smoothed out with a form factor (coupling function) 
g = g{k), k G M^, see (B.l) in Appendix B. The coupling constants Ai and A2 measure 
the strengths of the energy conserving (position-position) coupling, and the energy 
exchange (spin flip) coupling, respectively. Spin-flips are implemented by the Sj in 
(1.4), representing the Pauli matrix 
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(1.5) 



acting on the j-th factor of T^s ■ The total Hamiltonian takes the form 

H = Hs + Hn + v. (1.6) 

The dynamics of a density matrix pt of the system S-l-R is governed by the Liouville-von 
Neumann equation 

—pt = -i[H,pt], 

with initial condition pt\t=o = Po- The solution to the Liouville-von Neumann equation 
is given by pt = e^'*^poe'*'^- We are interested only in information on the subsystem 
S, so we trace out the degrees of freedom of R. The state of S is given by the reduced 
density matrix 

p, = TVR(e-«^poe"^), (1.7) 

where po is the initial density matrix of the coupled system, and Tr^ is the partial 
trace over the degrees of freedom of the reservoir. The operator acts on the Hilbert 
space T^s = C2 (g) . . . (g) C2 = C^"^ of S only. 
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Our goal is to analyze the time evolution of matrix elements of the reduced density 
matrix (1.7) in the energy basis, which plays a special role in quantum information 
theory. The energy basis consists of eigenvectors ip^ of H^, indexed by spin configura- 
tions 

a = {(Ji, . . . ,(TAr} e {+1,-1}''^, tpa = ^ai®---®^. 



Here, 



so that 





" 1 " 




' " 









1 



N N 

Hs^Pa = E{a)ipa with E{a) = ^ JijOiaj + BjOj 

i,j=l j=l 

We denote the reduced density matrix elements as 



(1.8) 

(1.9) 
(1.10) 
(1.11) 



The dynamics of the register alone (without coupling to the environment) is given 
by Pt = e~'*^SpQe'*^s^ where Pq = TrR,(po), so matrix elements of p^ have the time 
dependence 

We view the energy differences 



(1.12) 



TV N 

e{a,T) := E{a) - E{t) = ^ J,,(a,a, - r,r,) + ^i?,(a, - r,) (1.13) 

as being eigenvalues of the Liouville operator 

Ls = Hs^t-t^Hs, (1.14) 

acting on the doubled space 

Hs ® Hs = (C^ ® C^) • • • (g) (C^ ® C^), (1.15) 
where the j-th pair is the doubled space of the j-th qubit. 



Discussion of main results. In the resonance approach used in this work, 
we examine the influence of the interaction (1.4) on the free dynamics (1.12) for small 
coupling parameters Ai, A2. Under the perturbation, the phase factors e = E{T)—E(a) 
in (1.12) become complex resonance energies, Eg = ee{Xi,X2) £ C. The latter encode 
properties of irreversibility of the reduced dynamics of S (decay of observables and 
matrix elements - the dynamics of the entire system S + R is unitary, by contrast). 
We consider the regime where the resonance energies ee{Xi,X2) do not overlap as the 
perturbation is switched on, so that each resonance energy can be followed separately. 
This means that the coupling parameters must be small with respect to the gap between 
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differences of energies of Hs, see condition (Al) in Section 2 below. ^ Wc make as well a 
technical assumption (A2) on the regularity of form factors gi and g2 which we explain 
in Section 2. 

Dynamics of S. Our first result is a detailed description of the evolution of the 
reduced density matrix elements (and hence of all observables) . Set 

(([poo] a,r)) = lim ^ r[pt]a,r dt. (1.16) 
T-»oo 1 Jo 

We show in Theorem 2.1 that this limit exists, and that for all t >0, 

[PtU,J - (([Poo]_.,_r)) = E ' [ E * [PoU',^' + OiXl + Ai)" 

{e,s: 

+0((A? + Ai)e-'*), (1.17) 

where Imei*'' > and co' satisfies 2max{Imei'^^} < < r, with r > a constant 
depending on the regularity of gi, g2 (see Condition (A2) in Section 2, and also [16]). 
The * in the sum (1.17) means that we sum only over configurations a',T' such that 
e{_a',T') = —e. The coefRcents w are overlaps of resonance eigenstates (see Section 
2.1), which vanish unless e{a,T) = — e, in which case they are 0(1) in Ai, A2. The s'f^ 
are eigenvalues of a certain explicit operator K{ijj'), a "spectrally deformed Liouville 
operator" (see Section 3.1). They have the expansion 

e(f) = e + 5(:^+Oi\i + \i), (1.18) 

where the label s = 1, . . . , i^(e) indexes the splitting of the eigenvalue e of Ls, hav- 
ing multiplicity d{e), into z/(e) < d{e) distinct resonance energies. The lowest order 
corrections si^^ satisfy 

(5(^) = 0(A? + Ai). (1.19) 

They are the (complex) eigenvalues of an operator Ag, called the level shift operator 
associated to e (Ag is related to the Lindblad generator). This operator acts on the 
eigenspace of Ls associated to the eigenvalue e (a subspacc of the qubit register Hilbert 
space; see equation (3.21) for the formal definition of Ag). It governs the lowest order 

(s) 

shift of eigenvalues under perturbation. One can see by direct calculation that ImSi' > 
0.2 

Discussion of (1.17). To lowest order in the perturbation, the group of reduced 
density matrix elements [Ptl^.r associated to a fixed e = e{a,T) evolve in a coupled 
way, while groups of matrix elements associated to different e evolve independently. The 
density matrix elements of a given group mix and evolve in time according to the weight 

. (s) 

functions w and the exponentials e'*^« . In the absence of interaction (Ai = A2 = 0) all 

^Our method is applicable as well if this condition is not imposed. Work on this is in progress. 

^This can also be inferred from general considerations [14]: If the imaginary part was negative, then 
the average of some observables would explode as time increases, contradicting the fact that the total 
dynamics, a group of automorphisms, cannot increase indefinitely the average of any observable. 
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the Ee are real. As the interaction is switched on, the Eg typically migrate into the 
upper complex plane, but they may stay on the real line in certain cases. The matrix 
elements \fit]a,T of a group e approach their ergodic means (1.16) if and only if all the 

(s) 

nonzero Se have strictly positive imaginary part. In this case the convergence takes 
place on a time scale of the order 1/je, where 

7e = min|lm4^) : s = 1, . . . ,iy{e) s.t. e^^) / o| (1.20) 

(s) 

is the decay rate of the group associated to e. If an Eg stays real then the matrix 
elements of the coresponding group oscillate in time. A sufficient condition for decay 

(s) 

of the group associated to e is 7e > 0, i.e. Im^e > for all ,s, and Ai, A2 small. 

Decoherence rates. We illustrate our results on decoherence rates for a qubit 
register with Jij = (the general case is treated in Section 2.3). We consider generic 
magnetic fields defined as follows. For nj G {0, ±1, ±2}, j = 1, . . . , A'", we have 

N 

BjHj = <^ Uj = Vj. (1.21) 

Condition (1.21) is satisfied generically in the sense that only for very special choices 
of Bj does it not hold (one such special choice is Bj = constant). For instance, if the 
Bj are chosen independent, and uniformly random from an interval [i?mirn ^max] , then 
(1.21) is satisfied with probability one. We show in Theorem 2.3 that the decoherence 
rates (1.20) are given by 

7e = { + ^^^^^'^ + ' ^ n I + 0{\\ + A^). (1.22) 

t A2t/o, e = U J 

Here, yi is a contributions coming from the energy conserving interaction, uq and y2 
are due to the spin flip interaction. The term yi2 is due to both interactions and is 
of 0{\\ + A2). We give explicit expressions for yo; Vi-, 2/2 and 7/12 in equations (2.23), 
(2.13), (2.15) and (2.24). 

- Properties ofyi{e): yi{e) vanishes if either e is such that cq := Z^"=i(o"j — Tj) = 0, 
or the infra-red behaviour of the coupling function gi is too regular (in three 
dimensions gi cx \k\P with p > —1/2). Otherwise yi{e) > 0. Moreover, yi(e) is 
proportional to the temperature T. 

- Properties of y2{e): y2{e) > if g2{2Bj,T,) 7^ for all Bj (form factor 52 (^) = 
g2{\k\,'Fi) in spherical coordinates). For low temperatures T, y2{e) oc T, for high 
temperatures 1/2 (e) approaches a constant. 

- Properties o/yi2(e).' If either of Ai, A2 or eg vanish, or if gi is infra-red regular 
as mentioned above, then 2/12 (e) = 0. Otherwise yi2{&) > 0, in which case yi2(e) 
approaches constant values for both T — >^ 0, 00. 
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- Full decoherence: If 7e > for all e 7^ then all off-diagonal matrix elements 
approach their limiting values exponentially fast. In this case we say that full 
decoherence occurs. It follows from the above points that we have full decoherence 
if A2 7^ and g2{2Bj, E) 7^ for all j, and provided Ai, A2 are small enough (so 
that the remainder term in (1.22) is small). Note that if A2 = then matrix 
elements associated to energy differences e such that cq = will not decay on the 
time scale given by the second order in the perturbation (Af). 

We point out that generically, S + R will reach a joint equilibrium as t — 00, 
which means that the final reduced density matrix of S is its Gibbs state modulo 
a peturbation of the order of the interaction between S and R, see [16]. Hence 
generically, the density matrix of S does not become diagonal in the energy basis 
as t — > 00. 

- Properties of yo- Do depends on the energy exchange interaction only. This 
reflects the fact that for a purely energy conserving interaction, the populations 

are conserved [16, 17]. If g2(2Bj, S) 7^ for all j, then yo > (this is sometimes 
called the "Fermi Golden Rule Condition"). For small temperatures T, yo oc T, 
while yo approaches a finite limit as T — > 00. 

In terms of complexity analysis, it is important to discuss the dependence of on 
the register size N. 

- We see from (2.23) that yo is independent of A^. This means that the thermal- 
ization time, or relaxation time of the diagonal matrix elements (corresponding 

to e = 0), is 0(1) in N. 

- To determine the order of magnitude of the decay rates of the off-diagonal density 
matrix elements (corresponding to e 7^ 0) relative to the register size A^, we 
assume the magnetic field to have a certain distribution denoted by ( ) . It follows 
from the explicit expressions for yi, y2 and yi2 (see (2.13), (2.15) and (2.24)) that 



(yi) = yi oc el, (ys) = CB^ia - r), and (yi2) = cs(Ai, A2)Aro(e), (1.23) 



where Cb and cb = cb(Ai, A2) are positive constants (independent of A''), with 
cb(Ai, A2) = 0{Xi + A2). Here, No{e) is the number of indices j such that aj = tj 
for each {^,t) s.t. e{a,T) = e, and 



is the Hamming distance between the spin configurations a and t (which depends 
on e only). 

- Consider e 7^ 0. It follows from (1.22)-(1.24) that for purely energy conserving 
interactions (A2 = 0), 7e oc Xfe^ = Xl[J2f=i{<^j — Tj)]'^: which can be as large as 
0{XiN^). On the other hand, for purely energy exchanging interactions (Ai = 0), 
we have 7e oc A|lX) (_o- — r) , which cannot exceed ©(A^iV). If both interactions are 



N 




(1.24) 
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acting, then we have the additional term {yu), which is of order 0((Af + A|)A?^). 
This shows the following: 

The fastest decay rate of reduced off-diagonal density matrix elements due to the 
energy conserving interaction alone is of order \\N'^, while the fastest decay rate 
due to the energy exchange interaction alone is of the order . Moreover, the 
decay of the diagonal matrix elements is of oder Af, i.e., independent of N. 

- The same discussion is valid for the interacting register {Jij ^ 0), see Section 2.3. 

Remarks. 1. For A2 = the model can be solved explicitly [17], and one shows that 
the fastest decaying matrix elements have decay rate proportional to XfN^. Further- 
more, the model with a non-collective, energy- conserving interaction, where each qubit 
is coupled to an independent reservoir, can also be solved explicitly [17]. The fastest 
decay rate in this case is shown to be proportional to XfN. 

2. As mentioned at the beginning of this section, we take the coupling constants 
Ai, A2 so small that the resonances do not overlap (Condition (Al) in Section 2). 
Consequently \\N'^ and \'2N are bounded above by A = 2 minj=i^...,jv (see also 
Remark 4 after Condition (Al)) and thus the decay rates 7e do not increase indefinitely 
with increasing N in the regime considered here. Rather, the 7e are attenuated by small 
coupling constants for large N. They are of the order 7e ~ A. We have shown that 
modulo an overall, common (A^-dependent) prefactor, the decay rates originating from 
the energy conserving and exchanging interactions differ by a factor N . 

In this paper we prove the results only for sufficiently high temperatures. The 
general case will be treated elsewhere. 

3. The decay of off-diagonal matrix elements in the energy basis does not relate 
directly to measurements of entanglement, [18, 19]. We plan on elucidating the interplay 
between entanglement and decay of matrix elements in a subsequent work. 

Literature. Collective decoherence has been studied extensively in the literature. 
Among the many theoretical, numerical and experimental works we mention here only 
[1, 2, 4, 7, 8, 17, 20], which are closest to the present work. We are not aware of 
any prior work giving explicit decoherence rates of a register for not explicitly solvable 
models, and without making master equation technique approximations. 

2 Results 

As mentioned in the introduction, we assume that 

(Al) We have Co(|Ai| + jA2l)A^ < A for some constant Co (depending only on gi, g2, 
Jij and Bj). Here, A := min{e — e' : e, e' € S])ec{Ls),e ^ e'} is the gap in the 
spectrum of Lg. 

We implement a dynamical theory of resonances in a setting of spectral deformation 
(see Section 3.1). This leads to the following regularity requirement which we assume 
to be fulfilled throughout the paper. 
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(A2) The function form factors gi, 52 in (1-4) satisfy the fohowing condition. For 

h = gi 01 h = g2, 



is such that to hp{u + uj,a) has an analytic continuation, as a map C 
L^(M X S'^,du X da), into {|a;| < r}, for some r > 0. Here, (f) is an arbitrary fixed 
phase. 

Remarks. 1. Typically, the gap A depends on A^. We have ||-ffs|| < CN'^ and 
||Ls|| < CN"^, for some constant C. Therefore, if the 2^ (= dimT^g ® "^s) eigenvalues 
of Ls are roughly simple and equally distributed, then the gap A is of the order N^2~^ . 
In this case. Condition (Al) implies that the coupling constants Ai and A2 have to be 
exponentially small in the size N of the qubit register. However, the gap A tends to 
become larger as the multiplicities of the eigenvalues of Lg increase: A is the minimal 
distance between distinct eigenvalues of Lq, spread over an interval of size \\Ls\\. Due 
to the increase of multiplicities, the gap may become independent of A?^, as it happens 
in the following examples. 

- For Hs = J'^f=iSjSj_^_i (nearest neighbour interaction; and say Sf^_^-^ = Sf), 
we have spec(Ls) = J{-2N, -2N + 1, . . . , 2iV - 1, 2A^}. It follows that A = | J| is 
independent of A^. 

- For = X^jLi ^j^ji difference between two eigenvalues of Lg is given by 
e - e' = EiLi^iK - n'j)^ where nj,n'j G {-2,0,2}. Hence (for Bj > 0), A = 
2minj=i,...,iV-Bj. 

2. Examples of form factors satisfying (A2) are g{k) = hi{a)\k\Pe~^'^^ , where 
p = —1/2 + n, n = 0, 1, 2, . . ., and hi{a) = e^'^h\{a). They include the physically most 
important cases, see also [16, 17]. We point out that it is possible to weaken condition 
(A2) considerably, at the expense of a mathematically more involved treatment, as 
mentioned in [16]. The phase (j) has been introduced, and its physical interpretation 
has been given, in [9]. 

2.1 Effective dynamics of S 

The main result of this section is Theorem 2.1, in which we describe the effective 
dynamics of S and identify the dominant part. 

The evolution of reduced density matrix elements is governed by exponentials 
git4 ^{MM) ^ where £e''^(Ai, A2) are resonance energies, lying in the upper complex plane. 
The subindex e is the eigenvalue of Lg which the resonance branches out of: e^f' (0, 0) = 
e, and the index s = 1, . . . , i^(e) < d{e) distinguishes different resonance energies as- 
sociated to the same e {d{e) is the degeneracy of e as an eigenvalue of Lg). Using 
perturbation theory (we employ the Feshbach projection method (see Section 3.1 and 
[5, 16])), one obtains (1.18). 

Let {rli'^^}f-^]^ and {■qi'^^}^^]^ be bases of the eigenspaces of the level shift operator 
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Ag and its adjoint A* (see (3.21) for the formal definition of Ag), 

Aer?^'-) = ^r?^'-), r = l,...R, (2.1) 

^*~(s,r) ^ si%s,r)^ r = l,...R, (2.2) 

where R = R{e, s) is the geometric multiphcity of the eigenvalue di^^ of Ag. We choose 
bases that are dual to each other,^ meaning that 



'r/^'-WZ-'^'O = V- (2-3) 
We define the projection 

(2.4) 



Qe 

r=l 

acting on the eigenspace of Lg associated to e.^ 

Theorem 2.1 (Dynamics of matrix elements) Denote by /? the inverse tempera- 
ture o/R. There is a Xq > such that i/max{|Ai|, IA2I} < Ao//3 then the limit (1-16) 
exists for all a,T, and we have for t >0 

[Pt]a,r-{{[Poo]a,r))= (2.5) 



E 

{e,s: 4"' 5^0} 



+ R2{t). 



The * in the last sum indicates that we only sum over spin configurations o;',t' such 
that e{a' ,t') = — e. The remainders satisfy 

\Ri\ < CN^{Xj + Ai) and \R2(t)\ < CN^{Xj + Xl)e-'^'\ (2.6) 

where C is a constant, N is the register size, and where lo' satisfies 2 maxg 
lo' < t/2, with T given in Condition (A2). 

(s) 

Remarks. 1. Since is & projection with range in the eigenspace associated 

is) 

to the eigenvalue e of Lg, we have q\ Vt,ct = unless e{_a,T) = — e (see the scalar 
product in (2.5)). 

2. The condition max{|Ai|, IA2I} < Ao//3 stems from the particular complex defor- 
mation we choose in this work (translation). A mathematically more sophisticated 
treatment, involving a combination of spectral translation and dilation, and an itera- 
tive renormalization group analysis will yield the theorem for small Ai, A2, but with a 
temperature independent upper bound (see also [14, 15] and remarks in [16]). 

3. We mention again that in this work, we consider the regime of non-overlapping 
resonances, described by Condition (AI) at the beginning of Section (2). This means 
that Ai, A2 ~ l/N. 



^This is always possible, see Proposition A. 2 in Appendix A. 

''This projection is the same for all choices of bases 7)i"'^^ and rfe'^'' satisfying (2.3), as is easily 
verified using Proposition A. 2 of Appendix A. 



9 



2.2 Non-interacting qubit register in magnetic field 

We consider the qubit register Hamiltonian (1.1) with Jij = and Bj > 0, with a 
coupling to the reservoir given by (1.4). In this section we determine the resonance 
eigenvectors r]i'^'^\ rji'^'^ exphcitly, as weU as the resonance energies ee*-* to lowest order 
in the interaction, see Theorem 2.2. Those quantities are the key ingredients entering 
the dynamics which we describe in Theorem 2.4 below. 
Let a,T be spin configurations of the form (1.8). Then 

^a,T = VaiTi®---®^aNrN "^'^^^ (^q-t = <^(t ® V'r € (g) (2.7) 

is an eigenvector of Ls with eigenvalue e(_o",jr) = Y^jBjiuj — Tj). The genericness 
condition (1.21) implies that if V9o-,t and '-p a' ,r' a-re eigenvectors associated to the same 
eigenvalue, then Oj — Tj = a'- — rj for all j. If Oj — Tj = ±2 then Uj = ±1 and 
Tj = ^1 are determined uniquely, while if aj — Tj = 0, then there are two choices, 
aj = Tj = ±1. Consequently, an orthonormal basis of eigenvectors of Ls associated 
to a given eigenvalue e can be constructed as follows. Take any one eigenvector (pa,T 
associated to e and adjoin all linearly independent vectors a' ,t' with the property 
{(Tj — Tj = 0} \^a'j — Tj = 0}. Thus, with each eigenvalue e we associate the number 

No(e) = {number of indices j s.t. aj = Tj in any (^, jr) with e{_o;,T) = e}, (2.8) 

and the degeneracy of the eigenvalue e of Lg is d{e) = 2^°^^\ To each eigenvalue e of 
Ls there corresponds a unique sequence of -/Vo(e) indices indicating the locations j at 
which aj = Tj for all a, j associated with e. In other words, given e there is a unique 

sequence {lJ-k}k=i\ 

l<Hl<fl2<---<I^Noie)<N, (2.9) 

having the property that any eigenvector ^Pa,T associated to e satisfies 

jeifXk-. k = l,...,No{e)}. (2.10) 



Given an energy difference e (1.13), and a sequence g = {gj)f^i^\ Qj G {+1,-1}, 
we set 

Afu(e) 

Si^^ = A?[xi(e) + iyi(e)] + [x2{e) + iy2{e)] + (2-11) 



where 



xi(e) = -eoP.V.(5i,c^-V> Yl (^.12) 

{j: aj=Tj} 



ne, 



2 



yi(e) = ^7+, (2.13) 

X2{e) = - y aj F.y. [ u^g2i2u)coth{P\u\) — ^dn (2.14) 
{o- '^j¥=Tj} 

y2{e) = 217 Y B]g2{2Bj)coth{l3Bj), (2.15) 



3 



\ \^A^3 + 1) ± \j-h]{.Cj + \f + ^a\a-\hj{cj-\)\\ , (2.16) 
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with 



-A^eoP.V. (51, u; ^gi) , bj = AttX] 



(2.17) 



and 



eo 



JV 

= eo(e) = V((7,-r,), Qkiu) = \gk{\ul^)fd^, 7+ = lim t^gi(n). (2.18) 



The form factors gi, §2 (see (1.4)) are represented in spherical coordinates in (2.18) 
and P.V. stands for principal value. Note that eo is the same for all spin configurations 
a,T associated to the same energy e = e{a,T). This follows from the genericness of 
the magnetic field, (1.21), see paragraph after (2.7). We show in Theorem 3.5 that 



Imz^ > 0. Let us define the vectors 



ie) 
rie- 

Ve- 



AH 



where the ^af^.Tf^. at positions /ij, j = 1, 
given by 



■^ifar^r^, (2.19) 
■^^a^TN, (2-20) 
, iVo(e), are replaced by ^, ^ G (g) C^, 



?MjVo(e) 



it 



= f++ + 



1 + i 



1 + i- 



(2.21) 
(2.22) 



with normalization constant = [1 + bj ^c- ^{(bjCj — Imz^)^ + (a — Rez^)^}] ^. 

Theorem 2.2 (Resonance energies and states) Lete be an energy difference (1-13) 

and let Ag be the associated level shift operator. The vectors rji- and , (2.19) and 
(2.20), are bases of eigenvectors of Ae and its adjoint A*, respectively, which are dual 

to each other (see also (2.3)). The eigenvalues of Ke and A*^ associated to rji- andr^- 

are given by si- , (2.11), and its complex conjugate, respectively. Furthermore, we have 

=e + 5i^ +0{\\ + \^. 



Remark. The largest value of Nq is TV, which corresponds to e = 0, so d(0) = 2^. 
Here, jik = k = 1,...,N. The smallest value of A'^o is 0, which corresponds to 
e = zbcmax) where emax = 2 Bj is the largest eigenvalue of Lg. Thus emax is a simple 
eigenvalue of Ls- Here, no two aj, tj are equal, so the sequence {nk} is "empty". We 
have No{e) = No{—e), so d{e) = d{—e) for all eigenvalues e. 

The following result examines the resonance energies and shows expression (1.22) 
for the life times. 
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Theorem 2.3 (Fermi Golden Rule Condition and decoherence rates) Assume 
that the so-called Fermi Golden Rule condition is satisfied: 

Aiyo:=47rAi min {R2g2(25.) coth(/35.)} > 0. (2.23) 

j=l,...,N 

There is a c > s.t. if |Ai|, IA2I < c, then the decoherence rates are given by (1.22), 
with 

yi2(e) = ^ min|lmz+,Imz~| . (2.24) 
{j: aj=Tj} 

Remark. It is shown in Theorem 3.5 that Imz^ > provided abj 7^ 0, and that if 
a = 0, then z+ = 47riA|S|g2(2Gj) coth(^Bj), zj = and if bj = then zf = ±a. If 
abj = for all j then y 12(e) = 0. 

Let us illustrate how Theorems 2.1, 2.2 and 2.3 combine to give the detailed dy- 
namics of the register. Suppose that A2 / 0. It is clear that for generic values of the 
magnetic field, all si- are different for different g (see (2.11)). Thus all resonance en- 
ergies si- =6 + si- + 0(A^ + A2) are simple, for small enough Ai, A2. In this situation 
we obtain the following result: 



Theorem 2.4 (Dominant dynamics) Suppose A2 7^ and suppose that the mag- 
netic field is generic so that all SI-, (2.11), are distinct. There is a constant c s.t. if 
|Ai| + IA2I < c, then we have for all a, t 

[Pt]a,r-{{[Poo]a,r))= (2-25) 



where the * means that we sum only over spin configurations s.t. e{u' ,t') = —e, where 
ef^ = e + (5e~^ + 0{\\ + A2), the remainder terms Ri, R2 satisfy (2.6), and where 



Af„(e) 
j,k=l 

2.3 Interacting qubit register in magnetic field 

In this section we consider the Hamiltonian Ha,, (1-1), with generic parameters Jij and 
Bj. Energy differences of Hs are 

N N 

e{a, t) = E{a) - E{t) = i^i'^j ' ^^^j) + E ^3 - ^i) • (2-26) 

i,j=l j=l 
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0, 



(2.27) 



The condition e{a,T) = e{a',T') is equivalent to 

N N 
i,j=l j=l 

where rriij = a^aj — a[a'- — [riTj — t^Tj] and rij = aj — a'j — [tj — Tj]. For generic values 
of Jij and Bj, the only solution of (2.27) is m^- = 0, nj = for all i,j = l,..., N.^ 

Theorem 2.5 Let e be a nonzero eigenvalue of Lg. The resonance energies associ- 
ated to e are £e{o;,T) = e + de{^,T) + 0{Xf + A2), where {o;,t) varies over all spin 
configurations s.t. e{a,T) = e, and where (omitting {_£,t) in the notation) 

Se = Xfixi + iyi] + Xl[x2 + iy2] 

with xi{e) and j/i(e) given in (2.12) and (2.13), and 

N 



(2.28) 



1 



X2 



y2 



k=l 
N 



-^P.V. / u^\g2{u,T.)\ 



+ 



fc=l 



x52 



U + Vk 



Here, cq is given in (2.18) and 



Vk = -2a k 



N 



N 



'YiJjk + Jkj)Tj + Bk 



(2.29) 
(2.30) 

(2.31) 



The resonance eigenvectors associated to the resonance energy £e{^,T) are r?e(o-,T) = 

Va,T = Ve{a,r) (seC (2.1), (2.2)). 

This result shows that the decoherence rates induced by the energy conserving inter- 
action are again maximally 0(A^A^^), as in the case of the non-interacting register 
{Jij = 0). However, the decoherence rates induced by the exchange interaction have a 
complicated dependence on A'": y2 is a sum of A'" terms each one depending on N, the 
coupling parameters Jij and the magnetic field Bj. 

Remark. One can proceed as for the non-interacting register (Section 2.2) to 
analyze the resonances bifurcating out of the origin (determined to lowest nontrivial 
order by the spectrum of the level shift operator Aq). One finds that Aq has a simple 
eigenvalue at zero, and that the imaginary part of the smallest (nonzero) resonance is 
given by 



7o = 47rA9 min 

3=1,...,N 



q+g2(2C„+) ^ Cl_g2{2Cj,.) 



+ 



where 



N 

Cj,± = '^i'^jk + Jkj) =t Bj. 



(2.32) 



k=l 



^Indeed, to solve (2.27) with some rriij or rij nonzero means to introduce some correlations among 



the parameters Jij and Bj. Note that in particular, Jij 
parameters. 



J and Bj = B is not a generic choice of 
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3 Proofs 



3.1 Proof of Theorem 2.1 

In Theorem 3.2, we first obtain a suitable expression for the average {A)^ of an observ- 
able A G B{Hs)- This result is based on the dynamical resonance theory developed in 
[16], see also [10, 11], which we outline below. In a second step, wc carry out a refined 
analysis of the resonance theory to obtain Theorem 3.3. The combination of Theorems 

3.2 and 3.3 shows Theorem 2.1. 
Let AeB{ns). We have 

{A)t = TrsiPt A]=TTs+R[pt A0ln] 

= {ipo,e'*^^i'^2 [A(g)ls(S) Hr] e~^*^^i'^2,/;o) . (3.1) 

In the last step, we pass to the representation Hilbert space of the system (the GNS 
Hilbert space), where the initial density matrix is represented by the vector ipo (in 
particular, the Hilbert space of the small system becomes Hs <8) Hs)- 

The dynamics of an observable A is implemented by the group of automorphisms 
A ^ e'*'^^i'^2ylc^'*'^^i'^2 . The self-adjoint generator L\-^ \.^ is called the Liouville oper- 
ator. It is of the form i^Ai,A2 = -^0 + AiM^i + X2W2, where Lq = + Lr represents the 
uncoupled Liouville operator, and AiWi + A2W2 is the interaction (represented in the 
GNS Hilbert space). 

We take the initial state to be represented by the product vector ipQ = 7/^3,0 ® 
(the product form of the initial state is actually not necessary for our method to work, 
see [16]). Here, V's,o is an arbitrary initial state of S, and 'i/'R is the equilibrium state 
of R at a fixed inverse temperature < /3 < 00. We denote by V's,oo the trace state of 
S, (V'S,oo) (^S <8) lls)V'S,oo) = 2~^Tr (As). We introduce the reference vector 

Aef = -08,00 <X) tl^R. (3.2) 

The trace state has the separating property: given any state ips,o there is a (unique) 
operator B G SDTs, satisfying ^^3,0 = (Is ^ -B)V'S,oo- We write B' := ® B and note 
that B' commutes with all observables, so that we obtain from (3.1) 

{A)t = (^0, S'e'*^^i.^2 [A Is Ir] e-^*-^^i.^2V>,.ef> . (3.3) 

We now borrow a trick from the analysis of open systems far from equilibrium: one 
can find a (non-self-adjoint) generator i^Ai,A2 s.t. 

e"^^i.^2Ae"i*-f^Ai,A2 = e'*^^i'^2y4e""^^i>^2 for all observables ^4, t > 0, and 

i^Ai.AaV'ref = 0. 

There is a standard way of constructing Kx^^x^ given -Z^Ai,A2 a^id the reference vector 
ipref- ^\iM is of the form i^Ai,A2 = + Ai/i -|- A2/2, where the interaction terms 
appearing in the expression for Lx-^.x-z undergo a modification AiWi -|- X2W2 Ai/i -|- 
A2/2) c.f. [16]. Formally, we may replace the propagators in (3.3) by those involving 
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-^Ai,A2) ™d use that e ' ^i-^2ipj.^f = ^r^f- This procedure has been carried out in a 
rigorous manner in [16], yielding the fohowing resolvent representation 

{A)t = --^ / <^o, B' {Kx,,x,{u) - z)-' [A ® Is ® 1r] Vref) ^''dz, (3.4) 

where Kxj^^X2i'^) = Lq{uj) + Xili{uj) + A2/2('^)) Ii,2 are representing the interactions, 
and Lo 1-^ Kx^^Xii^) is a spectral deformation (translation) of Kai,A2- Relation (3.4) 
holds for < Imw < r (see condition (A2) in Section 2); the integrand is analytic in 
that domain, and continuous as Imw [ 0. For a; G M, the integrand is independent of 
CO and so it is constant for all lo in the domain of analyticity. 

The spectral deformation is constructed as follows. There is a deformation trans- 
formation U{(jj) = e~''^^, where D is the (explicit) self-adjoint generator of translations 
[16] transforming the operator -f^Ai,A2 ^ 

Kx,Mi^) = U{io)Kx,MU{io)-^ = Lo + u;N + Xih{Lo) + X2l2{u;). (3.5) 

Here, A'' is the total number operator of "Hr, having spectrum N U {0}, where is a 
simple eigenvalue (vacuum eigenvector iPr). For real values of lo, U{uj) is a group of 
unitaries. The spectrum of Kxi^Xii'^) depends on Imu and moves according to the 
value of Imcj, whence the name "spectral deformation". Even though U{uj) becomes 
unbounded for complex uj, the r.h.s. of (3.5) is a well defined closed operator on a 
dense domain, analytic in uj at zero. Analyticity is used in the derivation of (3.4) and 
this is where the analyticity condition (A2) of Section 2 comes into play. 

The point of the spectral deformation is that the (important part of the) spectrum 
of ArAi,A2('^) is much easier to analyze than that of Kxi,X2j because the deformation 
uncovers the resonances of A'ai,A2- We have 

spec {Ko{uj)) = {Ei - -Ej}i,j=i,...,Ar |J Wn + R}, 

n>l 

because Kq{uj) = Lq + ujN , Lq and N commute, and the eigenvectors of Lq = + 
are ipi ^Pj ipK- The continuous spectrum of Kq is bounded away from the isolated 
eigenvalues by a gap of size Imo;. The operator Ai/i(a;) -|- X2l2{oj) is infinitesimally 
small with respect to the number operator N, so for values of the coupling parameters 
Ai^2 small compared to Imuj, we can follow the displacements of the eigenvalues by 
using analytic perturbation theory. The following is an easy result (sec e.g. [14]). 

Theorem 3.1 Fix to' s.t. < u' < t (where r is as in Condition (A2) of Section 2). 
There is a constant cq > s.t. i/max{|Ai|, IA2I} < co/(3 ((5 is the inverse temperature) 
then, for all u with u' < < t, the spectrum 0/ A'ai,A2 ('^) ^he complex half-plane 
{Imz < w'/^} is independent of u and consists purely of the distinct eigenvalues 

{eW : e G spec(Ls), s = 1, . . . , v{e)], 

where 1 < u{e) < mult(e) counts the splitting of the eigenvalue e. Moreover, we have 
£e\^i-,^2) e as Ai,A2 0, for all e,s, and furthermore, Imei*^ > 0. Also, the 
continuous spectrum o/ii'Ai,A2('*^) ^^^^ region {Imz > 3a;' /4}. 
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Next we separate the contributions to the path integral in (3.4) coming from the 
singularities at the resonance energies and from the continuous spectrum. We deform 
the path of integration 2; = R — i into the line z = M + ia;'/2, thereby picking up the 

residues of poles of the integrand at e^^ (all e, s). Let ci^^ be a small circle around 
£e*\ not enclosing or touching any other spectrum of Kx^^x^i^)- We introduce the 
(generally non-orthogonal) Riesz spectral projections 

Q(^) = QW(a;, Ai, A2) = -^J^JK^.M^) - zy^dz. (3.6) 

It follows from (3.4) that 

u(e) 

= E E (^0' ^' ® Is ® iRlV-ref ) + R2, (3.7) 

e s=l 

where the remainder term R2 comes from the contour integral enclosing the continuous 
spectrum and satisfies 

\R2\ < CN^{Xl + A^)e-'^'*/2, (3.8) 

for some constant C not depending on the dimension N of Hs, nor on Ai, A2. Note 
that i?2 decays faster in time than each term in the main part. The estimate (3.8) is 
a direct consequence of Proposition 4.2 in [16] (see in particular equation (D.5) in the 
proof of this proposition). 

The ergodic mean time limits of R2 and all terms in (3.7) with ei^^ 7^ vanish, so 

{{A))^ := lirn^- J {A)tdt= ^ (V'o, 5' Q^'^fA ® Ir ® iRl^ref) • 

{e,s: ei'''>=Q} 

Combining the latter expression with (3.7) gives the following result. 

Theorem 3.2 For max{|Ai|, IA2I} < cq//? (see Theorem 3.1), we have for t > and 

{A)t-{{A))^= Yl ^''^^^\i^o,B'Qi'^[A0ls'S>MAe{) + R2, (3.9) 

{e,s: 

where R2 satisfies (3.8). 

Choosing A=\Lpr){ip^\ gives {A)^ = [pt]u,T- 

Theorem 3.3 Take max{|Ai|, IA2I} < cq/ (3 (see Theorem 3.1), let e he any eigenvalue 
of Ls and let a, t be .spin configurations. We have for t > 

E {fr',.',Qi'^^r,a) [Po]a',r' + 0(iV'(A? + A^), 

{^'iJt'- E(^')jr')=~^} 
where qi^^ is defined in (2.4), and N = dimHs is the register size. 
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We point out that Qe ipr,a vanishes unless e{T,a) = e. Theorem 2.1 now follows 
directly from Theorems 3.2 and 3.3. 

Proof of Theorem 3.3. Let R be the rank of qI"^ and [qI"^]*, and let {xe'''^}r=i 
and be bases of the ranges of those projections which are dual to each other, 

so that, by Proposition A.l, 

Qi'^=Elx^'^^)(x^'^^l■ (3-10) 

r=l 

We have 

i^M,A.(^)x^'')=4'^X^'') and [K^,,W^)Yx[''') = W^x^:^'\ (3.11) 

The following isospectrality result is inferred from the Feshbach method, sec [5], Section 
II and also [16]. We denote by the spectral projection of Kq associated to the 
eigenvalue e, and we set Pe = '^ — Pe- 

Lemma 3.4 (Feshbach map) Let x be an eigenvector of Kx^^x^{ijj) with eigenvalue 
£ (bifurcating out of e). Then ^ = PgX is an eigenvector of the operator 

Pe [e - I{uj)Pe(Kx„x,iuj) - e)-'PeI{u;)] Pe (3.12) 

with eigenvalue e. Conversely, if ^ RanPg is an eigenvector of the operator (3.12) 
with eigenvalue e, then 

X = [1 - Pe{Kx,M (^) - e)-^PeI{u)\ Pei (3.13) 

is an eigenvector of K\^^\^{u}) with eigenvalue e. Moreover, if x is an eigenvector as 
above, then ^ = PeX 7^ and conversely, if ^ is an eigenvector as above, then x given 
in (3.13) is nonzero. In particular, the geometric multiplicity of £ as an eigenvalue of 
^XiMi'^) is the same as that of s as an eigenvalue of (3.12). 

Expanding the resolvent in (3.12) around {Lq{u) — e)~^ we obtain for ^ = ^e*'*^"* 

^{^,r) ^ [^(.,r) ^ o{N\\l + V'R, (3.14) 

where rlf'"^'^ satisfies (2.1), with the level shift operator Ag defined in (3.21), and where 
N = dim?^s. We expand the resolvent in (3.13) around {Lq{ijo) — e)~^ and use (3.14) 
to obtain 

^{s,r) ^ _ p^(Xq(^) _ e)-iPe/(a;)] P, 4^'^^ Vr + 0{N\xI + Xl)) . (3.15) 

Proceeding in the same way we get the following representation for the eigenvectors 
Xe'^^ of the adjoint operator [Kxj^,\2{uj)]* , 

^ _ p^(L^(jj) _ e)-iPe(r )(a;)] Pe ^ Vr + 0{N\xI + Ai)), (3.16) 
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where r]i'^' satisfies (2.2). Relations (3.15) and (3.16) give 

=^(«,0^^j,+O(Ar(|Ai| + |A2|)), x^^) =^^'^)0^R+O(iV(|Ai| + |A2|)), (3.17) 

with the additional properties {1 ® (V'RlIXe*'^'* = fle'^'' ® + Oi^N'^{\\ + A^)) and 

{ICS (VrIIx^^) = t'^^ ® Vr + 0{N\\l + Ai)). 
Relations (3.10) and (3.17) show that 

(V'O, B' Q(^) [A® Is® iRjV'ref ) = (3.18) 
R 

J2 (V'O, B' 4''^^ V'r) {4''"^ ® V'R, [A® Is® iRlV'ref ) + 0{N\XI + AD) . 
r=l 

Let us take A = \ipT-){ip^\. Then we have (see also after (3.2)) 

(^(^''■), [A ® IslV's.oo) = 2^^/2 (^i^'^-), (^.,.) . (3.19) 

Next we insert a decomposition of identity, and use again the explicit form of ^ref) 
(3.2), to obtain 

^Vo, B' rii''"-^ ® Vr) = Yl <^0' ^' ® {'P^',r'^vi"'''^) 

= 2^/2 ^ (V;o,5' [|<^,,)(<^,,| lis ® iRlV'ref) (^a',r', ^ji^'"^^ ) 

= 2^/2 ^ (v:'.',r',r/(^''-)) [po]r',a'. (3.20) 

In the last step above, we commute B' to the right, B' jyi) {tp j-i | (8) Is (8) iRjV'ref ) = 
(V'o, [|¥'o-')(¥'r'| <8 Is ® IrIV'o) = [Po]t',ct'- Equations (3.18)-(3.20) demonstrate Theo- 
rem 3.3. This also concludes the proof of Theorem 2.1. □ 

3.2 Proof of Theorem 2.2 

The level shift operator associated to an eigenvalue e of Ls is defined as 

Ae = -PelPeiLo - 6 + iO)'^ PeIPe, (3.21) 

where is the spectral projection of Lq associated to {e}, Pe = t — Pg and Lq is the 
operator Lq restricted to the range of Pe-^ The interaction operator / has the form 

I = V- JA^/VJA^/^ (3.22) 

where J and A are the modular conjugation and the modular operator associated to 
the pair (9Jl, ^ref); where ^ref is the reference vector given in (3.2), and 9Ji is the von 
Neumann algebra of observables. The operator V is given by 

V = Xiivi (8) Is) (8 Mgi) + X2{v2 <8 Is) <8 M92)- (3.23) 
^Relative to [16], this definition differs by a sign. 
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Here, vi and V2 are the non-demolition and energy exchange interaction operators, see 
(1.4), and (p/^ig) = + ^tiio)) the thermal field operators. 

Theorem 3.5 (1) The level shift operator Ag has the form 

Ae = iA?yi(e) + Ai[x2(e)+iy2(e)]+ (^-24) 

{j: CTj=Tj} 

where yi{e), X2{e) and y2{e) are given in (2.13), (2.14) and (2.15), and where the 
operator is understood to act non-trivially only on the two dimensional suhspace, 

spanned by {ip^+,ip }, of the j-th factor (g in the Hilbert space (1.15). It is 

represented by the matrix 



= 

where a, b and c are given in (2.1 7) 



a + ibc —ibc 
—ib —a + 16 



(3.25) 

(2) The eigenvalues (2.16) of A'P satisfy lm.{z^) > 0, and they are strictly positive 
if ab 7^ 0. For a = 0, we have Zj' = ib{c + 1), = and for b = 0, we have z^ = ±a. 
The eigenvectors are given by (2.21). 



Since acts non-trivially on diff'erent factors of the Hilbert space for different j, we 
imm 
2.2. 



immediately see that the eigenvalues of Ag are the si-\ (2.11). This proves Theorem 



Proof of Theorem 3.5 The level shift operators (3.21), for interaction operators 
V = XG'Sits'SKpfsig) and reference states V'S./J^V'R (equilibrium state for the uncoupled 
dynamics), have been calculated explicitly in [16], Proposition 5.1. An easy adaptation 
to the present situation gives the following result. 

Proposition 3.6 We have the decomposition Ag = Ae,i -|- A|Ae_2; where A^ = 
limgj^o A#(e), with 

-2Ae,i(e) = Peivf^l-l®vf)Pe(gi,^^gi 

\ id + € 



-Pe{vi ® 11 -l^viYPe {^gi,cot\i{f3u/2)-^-^gij , (3.26) 

where lo > is the radial variable (spherical coordinates), and 
-2Ae,2(e) 

= Pe{v2 ® 1) / ""'If ^'I'f^ l-^S -e + u + ie)-\v2 ® t)Pe (3.27) 

+Peil^V2) [ ""'I'f ^'I'f r iLs-e + u + ie)-\l ^ V2)Pe (3.28) 
iMx52 |l-e+/^"| 

-Pe{v2^t) [ ""'I'f ^'""I'f r {Ls-e + u + ie)-\l V2)Pe (3.29) 

-Pe{l ® V2) [ ""'if ^'1_'g^?'' (i^S -e + u + ier\v2 ® l)Pe. (3.30) 
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All integration measures are dudS, where du is the Lebesgue measure, and dE is the 
uniform measure on S^. 

Remark. There are no "cross-terms" involving products of vi and V2 in the expression 
(3.21): Indeed, for instance 

AlA2Pe(^^l «) Is <8 (l>p{gi))Pe(Lo - e + iO)-^Pe{v2 (8) Is <8) M92))Pe = 0, 
since Sj (occurring in V2) maps any eigenspace of Lg into its orthogonal complement. 

Proposition 3.7 We have Ae,2 = aJ2(e) + iy2(e) + ire, where X2{e) and y2{e) are given 
in (2.14) and (2.15), and where 



{j: <Jj=Tj} 



with 



P = 47ri?2a2(25,) (l-Fj ^-^(^B,{Sim\ h _ ^-2i3B,{S'.m)\-\ 



(3.31) 



(3.32) 



The operator is understood to act non-trivially only on the two-dimensional subspace, 

spanned by (y? }, of the j-th factor (8) in the Hilbert space (1.15). The 

"flip operator" Fj is defined by Fjipa,T = Vain <8) • • • <8) V'(-a,)(-r,) <8) • • • <8) VajyTN- 
In the orthonormal basis {(/3++,<^ }, (see also (1-9)) the operator has the form 



-1 



_Q2(3Bj 
1 



(3.33) 



Proof. We leave out the tensor product symbols ® when no confusion should occur. 
Take a Va,T in the range of Pg. It follows from 



N 



that 

Ls('y2 <8) ll)(/Ca, 



N 

j=l \ ki^j 

N 

= ~ 2PjCrj) (fam ■ ■ ■ V{-aj)Tj ■ ■ ■ VaNTN- 



(3.34) 



We now apply Pe('y2 (8 1) = PeEfc=i('S'fc <8 1) to (3.34). The only contribution is 
coming from terms where = j in the resulting double sum: indeed, {S^Sj (8 l)(pa,T 
is orthogonal to the range of Pg unless k = j. It follows that 



(3.27) ip^,_ 



dudE 



u2|52(|u|,E)p ^ 



XS2 



\l-e-M 



^(-2P,a,+« + ie)-Va,. 



(3.35) 
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In particular, the operator (3.27) is diagonal in the energy basis. Proceeding in the 
same fashion one finds 

N 

Ls(l (8) V2)(Pa,T = + 2BjTj)ipa^ri ' ' ' fuji-Tj) ' ' ' V'untn (3-36) 

and 

(3.28) = / dndS ""'.'f ^'t'aS'' Y.^'^^n + « + ie)- Va,r. (3-37) 
JkxS2 |l-e+P"| ^ 

The operator (3.28) is thus also diagonal in the energy basis. Next we consider (3.29). 
We apply Pe{v2 ®t) = Pe Y.k=i{Sl 1) to (3.36). The only non-vanishing contribution 
comes from k = j va. the resulting double sum and only for terms where Gj — Tj = 0. 
We obtain 

«'b2(|n|,S)|2 

I UUUZj 

{j: aj=Tj} ' 

x{2Bjaj +U + fe)"Vairi • • • (p(-aj){-Tj) ■ ■ ■ H^ontn- (3-38) 
Note that aj can be replaced by Tj in (3.38). A similar argument gives 

«'b2(|«|,s)p 



(3.29) = - I d^dE^iMM^ 



(3.30) (^^,, = - J2 [ 



dudE 



c2 1 — e+/^" 

X {-2Bjaj + U + ie)~ Vain • • • (p(-aj){-Tj) ■ ■ ■ <faMTM ■ (3-39) 

The operators (3.29) and (3.30) are not diagonal in the energy basis. 

Next we use (3.27)-(3.30) and limg_^o+(-a + u + ie)~^ = -m5{u - a) + P.V. ^ 
in (3.35) and (3.37)- (3.39) to arrive at 

Ae,2<^a,T = [X2{e) +iy2{e)]ipa,T (3.40) 

E ^l n!e-2^i./ --.-- (3-41) 

We have |1 - e^^^^'^J]-^ = ^^^^5^^' so (3.41) plus (3.42) combine to 

4i7r B]g2{2B,) {1 - F.e-'^^^'^^} _ Jl^pB.a.f ^,!- (3-43) 

{j: aj=Tj} 

The form (3.33) of in the basis {(p^^,ip } is immediately obtained from (3.32). 

This completes the proof of Proposition 3.7. □ 
The following result follows directly from (3.26). 
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Proposition 3.8 We have 

Ae,i = 12/1 (e) -eoP.V. (51,0;- V> Yl ^^-^^^ 

{3: (Jj=Tj} 

where yi and eo are given in (2.13) and (2.18). 

We obtain (3.24) now by combining Propositions 3.7 and 3.8. This shows point (1) of 
Theorem 3.5. Point (2) is verified easily by using the expression (2.16). □ 

3.3 Proof of Theorem 2.3 

According to (1.20) we have 7e = min{lmei-^ : q G {+1,-1}^o(^) s.t. ei" / }, 
where e = e(a, r) = Y.^^^ Bj{aj - tj). 

- For e 7^ and Ai, A2 satisfying (Al), we have e + si- ^ 0, and hence ei- 7^ 0, 
if |Ai|, IA2I < c, for some c > 0. The smallest imaginary part of si- (for e fixed) is 
A?yi(e) + Aiy2(e) +yi2(e). 

- For e = 0, we have e + di- = <Tj=Tj} ^j^ ■ Indeed, e = forces aj = tj for all 
j, and eo = (see (2.18)). It follows that aj = and so 

zf = ATTi\lB]g2{2Bj)coth{PBj) and z-=0. (3.45) 

The smallest imaginary part of di- is thus zero, corresponding to gj = —1 for all 
j = 1, . . . , Nq{Q) = N.''' All other imaginary parts are strictly larger than the gap given 
by (2.23). 

This shows formula (1.22) and completes the proof of Theorem 2.3. □ 

3.4 Proof of Theorem 2.5 

Let (pa,T be an eigenvector of Ls associated to the eigenvalue e{a,T). Let k = 1, . . . ,N 
be a fixed index. The vector {S^<SiS^)ipa,T is again an eigenvector of Lg with eigenvalue 
e(o-', r'), where {crj,rj) = {crj,Tj) for all j 7^ k, and (o-j^,T^) = (-o-fe,-Tfc). We now 
show that e{_a,T) 7^ e{a',T') unless e{a,T) = 0. Indeed, suppose that e{_a,T) = 
e{_a',T'). Then, due to the genericness of the parameters Jij and Bj (see after (2.27)), 
we have = 0, from which it follows that Ufc = r^. Furthermore, since niik = for all 
i, we obtain ai = Ti for i = 1, . . . , N . We conclude that for all k, 

n(a,r)(5^®5fcVa,r = 0, ife(a,T)7^0, (3.46) 

where Pg is the spectral projection of Ls associated to e. One sees also easily that if 

k^l, then Peia,r){S^k ® Sf)cp^,r = 0. 

^It can be inferred from general considerations that at least one eigenvalue of Ao must be zero. 
Indeed, since the generator i^Ai,A2 has been designed to annihilate the reference state iprei = i^s,oo'S>ipR, 
it follows that AoV's.oo = [12]. Note that indeed, for Qj = —1, all j, we have tJq-' = V's.oo- 
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Proposition 3.9 The level shift operators Ag with e 7^ are diagonal in the energy 
basis. Their eigenvalues are given by Xi[xi +iyi] + X2[x2 + 12/2] (see (2.12), (2.13) and 
(2.29), (2.30)). 

Proof. The spectrum of Ae,2 (see Proposition 3.6) is obtained as in the proof of 
Proposition 3.7. Relations (3.34) and (3.36) are replaced by 

N 

Ls{v2®t)ip^^r = ^[e{a,T)+Vj{a,T)\ip„^ri---'f{-ai)rj---'f<TNrN (3-47) 
N 

where Vj, Vj are given in (2.30). The terms (3.38) and (3.39) vanish due to (3.46). It 
then follows easily that the spectrum of Ae,2 is [x2 + 1^2] • 

The operator Ag^i is the same as in the case Jjy = 0, so we can again use Proposition 
3.8, and, together with Proposition 3.6, this gives the result. □ 

A Dual bases, projections, resonance eigenvectors 

Proposition A.l Let Q be a finite- dimensional projection in a Hilbert space H. Given 

any basis {Xr} ofHanQ, there is a unique basis {Xr} of Han Q* satisfying the duality 
condition {Xr,Xr') = S^y . For Xr> Xr obtained in this way, we have Q = |Xr)(Xr|- 

Proof. Take any basis {xr} of RanQ and let € TC be arbitrary. We have Qip = 
'^r^rCri'ip), whcrc ip Cr{ip) G C is a linear functional. Consequently, for each r, 
there is a £ such that Cr(V') = {Xr, V')- Hence Q = \Xr){Xr\- The vector Xr' is 
left invariant by Q, so it follows that {Xr,Xr') = ^r,r'- 

We show that the Xr are a basis of RanQ*. Firstly, wc have Q*Xr = Xr since 
Q* = X^,. |Xr)(Xr| and {Xr,Xr') = ^r.r'j SO we Only need to show linear independence. 
Let Zr be scalars. If ZrXr = 0, then, by taking the inner product with Xr', where r' 
is arbitrary, and using that {xr, Xr') = ^r,r', we see that z^/ = 0. 

To show uniqueness of {x^} foi' fixed {Xr}) we suppose that {5, } is another dual 
basis. Then aj. = fJ,r,r'Xr' and by the duality condition, (Xr^Q^r') — ^r,r' — 

J SO 

ji is the identity. □ 
Proposition A. 2 There are bases {r]i'^^}r o^nd {'ffi'^^}r of the eigenspaces of and 

is) 

(Ag)* associated to the eigenvalue 5e and its com.plex conjugate, satisfying the du- 
ality property (2.3). Those bases are unique in the following sense: any other pair 
of such bases {oe^''^},., {Se*'^''}^. is given by ai^'^"^ = Ylr'\.-^]r,r''ne"^ ^ (^f^d a'f'^^ = 
Ylr'[(^~^T]r,r''ne''' \ whcrc A is an invertible matrix. 

Proof. Prom Proposition A.l we know that we can find bases {xe^'^^}r and {xe'^^}r 
of the eigenspaces of Kxj^,\2{uj) and its adjoint, associated to the eigenvalue ei^^ and 
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its complex conjugate, respectively, so that \x'e'^\xe'^ / = ^rr'- Expansions (3.15), 



(3.16) show that (rii'''\7f/'''^) = \imx„x,-.o {xi'"\xi'"'^ ) = Or,.' 

We know from the proof of Proposition A.l that there is a unique dual basis for fixed 
{il'i'''^}. Thus, any pair of dual basis is gotten by a change of basis of one particular 
such pair. Let ai'^^ be obtained by a base change matrix A as in the proposition. 
There is a unique associated dual basis 5^*'*^ = X^^/ [-Bj^^^'^e*'^ ^ ■ It is easy to see that 
ai'''\ai'''''^) = 5ry implies that B = {A-^)*. □ 



B Operators Kx.^x^ and Kx^,\^{uj) 

The purpose of this Appendix is to provide some details on explicit formulas of the 
operators i^Ai,A2 Kx^^x2{^)- For more detail, we refer to [16]. 
Smoothed-out creation and annihilation operators are defined by 

a*ig)= [ g{k)a*{k)d% a{g) = f ~^a{k)d?k, 

for g = g(k) G L^(M^,d^/c), and the field operator is given by 

^{g) = l=[a*{g)+a{g)]. (B.l) 

The so-called Araki- Woods representation gives the Hilbert space (GNS) representation 
of the infinitely extended Bose gas in thermal equilibrium [3, 13].® The Hilbert space 
is given by the bosonic Fock space over the one-particle space L^(M x S'^,d?k x dS), 

jr = jr(L2(Mx52,d2A;xdS)). (B.2) 

The thermal annihilation operators are 



0/3(/) = + I^l3{u)x+{u)uf{u,a)) - a*{e"i'<^ fifs{-u)x-{u)uf{-u,a)), (B.3) 

where iJ,p{u) = (e''" — 1)~^, x± ^re the indicator functions of R±, and G M is an 
arbitrary phase. The a*p{f) are obtained by taking the adjoint on the r.h.s. of (B.3). 
It is easy to see that the OCR are satisfied. The thermal field operator (B.l) is thus 
represented by 

<j>pU) = -^{a}{f) + ap{f)) = -L(a*(/^) +a(/^)) =: ^(/^), (B.4) 

for / G L^(R^), where is defined in (A2), and where the cp iii the r.h.s. is the 
field operator in JT. The equilibrium state is represented by the vacuum vector of 



*In this paper, we directly work in a spatially unitarily equivalent representation of the original 
representation, see [16] for details. 
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^R,/3 = ^- For a one-body operator O acting on wave functions of the variables {u, a), 
we write 

dr(0) = / a*{u,a)Oa{u,a)duda. (B.5) 

for the second quantization of O. The dynamics of the field is generated by 

Lr = dr(u), (B.6) 

the second quantization of the operator of multiplication by u. We have L-^Q^^p = 
0, and for z G C, e^^«0^(/)e-^'^"» = 2-^2 [a^^e"^"/) +a^(e^"/)], which gives the 
dynamics for z = it. 

The Liouville operator -Z^Ai,A2 acting on (C^ (gi C^) (g) ^ is given by 

Lx,,x, = Lo + XiWi + X2W2, (B.7) 
Lo = Ls + LK = ifs®]ls-]ls®^fs + dr(n), (B.8) 

AT 

^fe = J2^J^'^S^'f''^^9k)p), k = l,2, (B.9) 

i=i 

where we understand Sj = Sj and Sj = Sj. 

The deformation group U{io) (see after (3.4)) is the translation group U{io) = 
g-ia;dr(ia„)^ and the spectrally deformed Liouville operator is 

Lx,,x,{uj) = Lo+uN + XiW{uj) + X2W2{uj), (B.IO) 

where N = dT(l) is the number operator in JF, and where Wk{oo) = e^'^'^^^^^'^Wke'^'^^^^''^ 
(see also (B.18)). 

Definition of the operator K\-^^\^. This operator can be expressed in terms of the 
non-interacting Liouville operator Lq, the interaction AiWi + X2W2, sec (B.7)-(B.9), 
and the modular data J, A associated to the vector Vref) (3-2), and the von Neumann 
algebra 9Ji = B{Tls) Is ® where 9Jl;3 is the Weyl algebra of the Bose field (see e.g. 
[6, 16]). J is an anti-unitary operator and A is a self-adjoint non-negative operator. 
The defining properties of J and A arc JA^^'^Mftg^Q = M*Jlfl 0, for any M G 9Jt, where 
M* is the adjoint operator of M. The explicit expressions are (see also [6, 16, 15]) 

J = Js0Jr and A = As Ar, (B.ll) 
As = e-f^^^, (B.12) 
Ar = e-^-^i^, (B.13) 

Js(l)l(g><pr = C<pr(^C(pl, (B.14) 

jRipniui,cri, . . . ,Un,an) = e''^'^ip^^{-ui, ai, . . . , -Un, an) , (B.15) 

where the action of the antilinear operator C is to take the complex conjugate of vector 
coordinates in the basis {ipj}jLi of Hs, and ■i/'n is the complex conjugate of ipn G J^. 
Relation (B.15) shows that JRa*{f{u,a))jR = a*{e''t'J{-u,a)), for / G L'^{R x S^). 
The interaction operators Ik in (3.5) are given by Ik = Wk — Wj^, where 

^ 1 r 1 

W'k = JA'/'WkJA'/^ = Is ® ^ 5| ® -= [a* {{gk)p) + a{e-^-{gk)(3)\ • (B.16) 

.7=1 ^ 
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The spectrally deformed operator Kx^^^x^^uj) is obtained as follows. The transformation 
of creation and annihilation operators under U{u), (B.IO), is 

U{u)a*{f)U{u)-^ = a*{f{- + cj)), a; G R, (B.17) 

where /(• + oo) is the shifted function {u,a) i-^ f{u + ui,a). Relation (B.17) can be 
written in the form U {uj)a!^ {f)U = a*(e'^^"/). In order to obtain an analytic 
extension of (B.17) to complex lo, we need to take the complex conjugate of uj in the 
argument of the annihilation operator (since the latter is anti- linear in its argument). 
We thus have /fc(w) = Wk{uj) — W^(w), with 

N 



Wk{i^) = Y,S^^®ts®^[a\{9k)p{- + oj))+a{{gk)p{- + uJ)% (B.18) 
j=l V 

^ 1 r -1 

Wi{u) = ]ls®^5|®-= a*{{gj,)0{- + u)) + a{e-^^^+'^\gk)0{-+uJ)) . (B.19) 
i=i V2 



3 

Finally, we have Kx^^X2{^) = ^x^Xii^) + Ai/i(w) + \2I2iu}). 
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